Scale-Locality of Magnetohydrodynamic Turbulence 



o 
o 

(N 
O 

Q 

00 



43 : 

Oh- 
i < 

o 

CO 

o3 



> 

in 

m 

On 
O 



13 



Hussein Aluie 1 ' 2 and Gregory L. Eyink 1 
1 The Johns Hopkins University, Applied Mathematics & Statistics, Baltimore, MD 21218, USA 
2 Theoretical Division (T-5/CNLS), Los Alamos National Lab, Los Alamos, NM 87545, USA 

(Dated: December 18, 2009) 

We investigate the scale-locality of cascades of conserved invariants at high kinetic and magnetic 
Reynolds numbers in the "inertial-inductive range" of magnetohydrodynamic (MHD) turbulence, 
where velocity and magnetic field increments exhibit suitable power-law scaling. We prove that fluxes 
of total energy and cross-helicity — or, equivalently, fluxes of Elsasser energies — are dominated by 
the contributions of local triads. Corresponding spectral transfers are also scale-local when defined 
using octave wavenumber bands. Flux and transfer of magnetic helicity may be dominated by non- 
local triads. The magnetic stretching term also may be dominated by non-local triads but we prove 
that it can convert energy only between velocity and magnetic modes at comparable scales. We 
explain the disagreement with numerical studies that have claimed conversion nonlocally between 
disparate scales. We present supporting data from a 1024 3 simulation of forced MHD turbulence. 

PACS numbers: 95.30.Qd, 52.35.Ra, 47.27.Jv 



Magnetohydrodynamic (MHD) turbulence is pervasive 
in astrophysical systems. Turbulent plasma fluctuations 
commonly possess power-law spectra over vast ranges 
of scales where both viscosity and resistivity are neg- 
ligible. We call such ranges "inertial-inductive", since 
nonlinear dynamics (inertia/Lorentz force and convec- 
tion/induction) dominates the physics at these scales. 
MHD plasma turbulence, with power-law scaling of both 
spectra and structure-functions in the inertial-inductive 
range, plays a central role in star formation, accretion 
of matter near active galactic nuclei, solar physics, and 
in the generation of large-scale magnetic fields in such 
systems. There are several competing theories for the 
spectrum of strong MHD turbulence, including those of 
Iroshnikov-Kraichnan [l|, Q, Goldreich-Sridhar and 
Boldyrev [4( . All of these theories assume scale- locality 
of the nonlinear cascade, following the classical ideas of 
Richardson, Kolmogorov and Onsager for turbulence in 
neutral fluids. Scale-locality is fundamental to justify the 
universality of the postulated turbulent scaling laws. 

A consensus has been forming in recent years, however, 
that cascades in MHD turbulence are nonlocal processes 
[H> S HI • Schekochihin et al. [B[ emphasized the non- 
local nature of the interactions between the velocity and 
magnetic fields as a hallmark of isotropic MHD turbu- 
lence. This conclusion was reaffirmed in several subse- 
quent studies, most categorically by Yousef et al. 
who claimed that there is a direct exchange of energy 
between motions at the largest scales in the system, at 
which the flow is being stirred, and the magnetic field at 
arbitrarily small scales in the inductive range. A more 
refined analysis of the locality of interactions was carried 
out by Alexakis et al. Q who concluded, based on direct 
numerical simulations (DNS) of MHD turbulence that 
the magnetic field gains energy at scales £ in the inertial 
range from the straining motions at all larger scales > £ 
and especially from the forcing scale L>1 Carati et al. 
subsequently carried out DNS at a higher resolution, 



arriving at conclusions similar to and, furthermore, 
claiming that there is non-local transfer of Elsasser ener- 
gies as well. Related ideas have surfaced in the accretion 
disk community [^, [13] • 

In this letter we address the scale-locality of MHD cas- 
cades by a direct analytical study of the equations: 



d t n + (u ■ V)u = -V P + (b • V)b + ^V 2 u + f 



d t b + (u • V)b = (b-V)u + r/V 2 b 



(1) 



for V ■ u = V ■ b = 0. Here b = B/^/47rp is the magnetic 
field in Alfven velocity units and p is total pressure 
(including magnetic pressure). Our main conclusion 
is that, under very weak scaling assumptions, MHD 
turbulence has scale-locality properties only a little less 
robust than those of hydrodynamic turbulence. We will 
support our analysis with a pseudospectral DNS at 1024 3 
resolution with phase-shift dealiasing. For our numerical 
work we choose viscosity v and resistivity r\ to be both 
equal to 1.1 x 10~ 4 . The external stirring force is a 
Taylor-Green flow f = /o[sin(fc/:r) cos(fc/y) cos(fc/z)x — 
cos(fc/a;) sin(fc/y) cos(fc/z)y] applied at modes 
kf = 2 with an amplitude /o = 0.25. The 
Reynold's number based on the Taylor scale X u = 
2ttV£^/( / dkk 2 E u (k)) 1/2 is i?e A „ = u rms \ u jv = 909. 

Our proof of local cascade of invariants in MHD 
turbulence is very similar to that given for hydrody- 
namic turbulence in 
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Il2l 1 1 31 ] . We employ the spa- 
tial coarse-graining approach, commonly used as a mod- 
elling tool in the Large-Eddy Simulation (LES) com- 
munity fl5j |. Coarse-grained fields are defined by 
ft(x) — J dr G^(r)/(x + r), with a filtering kernel 
Gi(r) — £~ 3 G(r/£) which is sufficientl y s mooth and 
decays sufficiently rapidly for large r [11(. Coarse- 
grained MHD equations can then be written to de- 
scribe U£ and hi, along with corresponding budgets 
for the quadratic invariants — energy, cross-helicity, and 
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magnetic- helicity — at scales > £. See [16| . E.g. the 
time-derivative of large-scale energy (l/2)[|u^| 2 + |t>£| 2 ], 
in addition to space-transport terms, contains also as sink 
terms the kinetic energy flux —II" = Vuf : Tg and the 
magnetic energy flux —II* =j e -eg, withj £ = Vxbf. Here 
T£ t ij — Ti(ui,Uj) — Ti(bi,bj) is the total stress generated 
by scales < £, both the Reynolds stress and the Maxwell 
stress, and se,i = CijkTe(uj, bk) is the electromotive force 
generated by scales < £. We employ the notation 



n(f,g) = (fg)i - f e ge 



(2) 



for the "central moments" of any fields /(x),<?(x) [14j . 

There are two facts crucial for scale-locality of the en- 
ergy fluxes Hg' b . First, all the filtered gradient-fields and 
the central moments can be expressed in terms of incre- 
ments. In general, for any fields, 



Vf e *5f (£)/£, n(f,g)ttSf(l)Sg(t), ti~-5f(£) 



(3) 



where increments are <5/(x, r) = /(x + r) — /(x), Sf(£) = 
sup r< ^ |<y/(r)|, and f' t = f — f e is the fine-scale (high- 
pass filtered) field. For details, see [HI]. The second 
crucial ingredient for locality is the scaling properties of 
the increments of velocity and magnetic field: 

Su(£) ~ f Sb(£) -to, < a u<b < 1, (4) 

where these relations may be assumed to hold either 
pointwise, with a the local Holder exponent, or in the 
sense of pth-order means, |j(5/|| P = (\8f(£)\ p ) 1 ^ p ', with a 
equal to l/p times the scaling-exponent £ p of the pth- 
order structure function. As long as < a u .b < 1, then 
(either locally or in the i p -mean sense) the fluxes n"' b 
are determined by modes all at scales comparable to £ 
[ill ]. F° r example, the contribution to any increment 
8f(£) from scales A > £ is represented by Sf A (£). Since 
the low-pass filtered field f A is smooth, its increment 
may be estimated by Taylor expansion and ((3]),(j4|) as 



6f A (£)~HVf A )~£A° 



r(£/A) ] 



and this is negligible for A ^> £ as long as a < 1. On the 
other hand, the contribution to any increment Sf(£) from 
scales S < £ is represented by 8f' s {£). Since f' s ~(Sf)(S) 
(even without taking any difference), (g]) implies that 



Sf's(i) 



£ a (6/£V, 



and this is negligible for i5 <C £ as long as a > 0. 

It is important to emphasize that the scaling laws like 
used in our proof are obtained in all theories of strong 
MHD turbulence. The Iroshnikov-Kraichnan theory pre- 
dicts that a u = (ib = 1/4. The Goldreich-Sridhar the- 
ory predicts distinct scaling for increments with displace- 
ments in different directions relative to a background field 

1/2 

b , with 5u(£\\) ~ &b(£\\) ~ £,/ for displacements in the 



(A/3 



for 



field-parallel direction and Su(£± 
displacements in the perpendicular direction. Such dis- 
tinctions make no difference to our proof, so long as both 
exponents an, a± lie between and 1. Similarly, our 
proof is fully compatible with possible intermittency cor- 
rections to scaling exponents. Although the precise scal- 
ing of stron g M HD turbulence is an open issue, numerical 
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sup- 



simulations , |18| and natural observations 
port the validity of the weak condition (j4|) for sufficiently 
high kinetic and magnetic Reynolds numbers. 

Our arguments imply also the scale-locality of cascades 
of the Elsasser energies (l/2)|z ± | 2 , with z ± = u =p b. 
This may be seen by considering the time-derivative of 
the large-scale energy densities (l/2)|z^| 2 , for which the 
sink terms are the fluxes 



V^:^(zT )Z ±) 



Sz T (£)[Sz ± (£)] 2 /£. Since these fluxes are expressed in 
terms of increments, they are scale-local under the weak 
condition (j4|). This may also be seen from an alter- 
native expression for the Elsasser energy fluxes which 



_ follow from the Politano-Pouquet relations [2l|, n 



-$/4£)(t'8x?(t)\5*F(t)Y 



where 



/ ang 



denotes av- 



erage over the displacement directions I. The scale- 
locality of cascades of the Elsasser energies is particu- 
larly important since the foremost phenomenologies of 
strong MHD turbulence P, [3, d, Q are based on the pic- 
ture of counterpropagating Alfven wavepackets expressed 
by the Elsasser variables . In terms of these vari- 
ables, the scale-locality properties of MHD turbulence 
are essentially the same as those of hydrodynamic tur- 
bulence. Scale-locality of the cascades of Elsasser en- 
ergies implies scale-locality of the flux of cross-helicity 
Tie-be — (l/4)|z/| 2 — (1/4) |z £ ~ | 2 (as well as scale-locality 
of flux of total energy (l/4)|z + | 2 + (l/4)|z^| 2 ). 

One cascade in MHD turbulence which may be essen- 
tially different is that of magnetic helicity. The time- 
derivative of large-scale helicity-density hi-He (where 
a> = (curl) -1 ^), in addition to space-transport terms, 
contains as a sink term the magnetic helicity flux — H 1 ^ = 
2bf£(.. Although £( ~ Su(£)5b(£), the coarse-grained 
magnetic field he will generally be dominated by modes 
at the forcing scale L. Thus, magnetic-helicity flux may 
possibly be dominated by non-local triads, with one mode 
at the large scale L. Similar issues arise for magnetic line- 
stretching. The time-derivative of large-scale kinetic en- 
ergy (1/2) | | 2 , in addition to space-transport terms and 

the sink term — n™, contains —h e Sfhg where the matrix 
Se — (l/2)[( VHf) + (VlIf) T ] is the strain from scales 
> £. Likewise, the time-derivative of large-scale mag- 
netic energy (1/2) | | 2 , in addition to space-transport 



terms and the sink term — In, contains 



-b^ S e bi 



Thus, 



this term represents conversion between large-scale ki- 
netic and magnetic energy by stretching of coarse-grained 
field-lines. Just as for magnetic helicity flux, this is a "hy- 
brid" quantity with both energy-range and inertial-range 
components. Although ~ 6u(£)/£, the coarse-grained 
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magnetic field hi can be dominated by modes at scale L. 
Thus, we cannot conclude that this quantity is dominated 
by local triads with all modes at scale £. 

Indeed, much of the recent discussion about apparent 
non-locality in MHD turbulence has revolved about this 
conversion term. One of the startling claims that has 
been made in recent numerical studies [1, 0, @| is that 
conversion between kinetic and magnetic energies pro- 
ceeds very non-locally, with magnetic modes at scale £ 
gaining energy equally from all velocity modes at scales 
> £ or even predominately from scale L ^> £■ In order 
to examine this claim, we must refine our methodology 
to consider band-pass energies. Following [12 ]. we define 
pointwise the kinetic and magnetic energy densities in 
the interval of scales \£, £] for £ > £, as 



e^, = (l/2)r^,H), efo = (l/2) 



r[Di,bi) 



Note that (•) now denotes scale £ and (•) scale £. Their 
time-derivatives are easily calculated to be 



d t e 



-(bibjSi 



bibjSij) + 



hbjSij / 



n - (rr) 



n 



(n b ; 



(5) 



(6) 



where • • • denotes total divergence terms that correspond 
to space-transport. As before, II" = — S:t and II = 

— J-e, and note that the double-filtering length-scale £ w 
£ for £ ^> I. It is "obvious" from these equations that 
the magnetic stretching terms transfer energy between 
velocity and magnetic-field modes only within the same 
band of length-scales \£, £} . Clearly, whatever energy is 
lost or gained from one field by line-stretching reappears 
in or disappears from the other field at the same scale. 
Non-colliding Alfven waves are an example of such non- 
local triadic exchange which is mediated by a uniform 
magnetic field at the largest scales, but which does not 
contribute to energy transfer across scales. 

Our conclusion above requires some caution, however. 
A counterexample is the Batchelor (viscous-inductive) 
range that occurs in MHD turbulence with a large mag- 
netic Prandtl number Pr m = v/r] ^> 1 [23]. This range 
consists of length-scales £ v 3> £ ^S> £ n far below the 
inertial- inductive range L 3> £ 2> £ v , with £ v and £ v the 
viscous and resistive length-scales, resp. In the Batchelor 
range, the energy is transferred directly from the veloc- 
ity modes at the viscous scale £ v into the magnetic-field 
modes at scales £ <C £ v . To see that this follows from 
our eqs. ©-([6]), we observe that the velocity-gradient 
in the Batchelor range is almost spatially constant and 
Vui « Vu for all It is thus easy to see that the 

stretching term in ([5]) equals —SijT(bi,bj) whereas the 

two flux terms become Sy (t(6,, bj)—r(bi, bj)). (Note that 



the stress in this range is almost entirely Maxwellian). 
These terms exactly cancel, by the Germano identity 
14, [ill. Thus, the line-stretching term acts as an ef- 



fective source to magnetic energy e - ~ , supplied by the 

flux of kinetic energy directly from the viscous scale £ v . 

The moral of this example is that the energy fluxes 
also contain line-stretching effects which must be consid- 
ered. Nevertheless, our conclusion is not altered that, 
in an inertial- inductive range, energy conversion by line- 
stretching is between velocity and magnetic-field modes 
at similar scales. The key point here is the scale-locality 
of the fluxes, which has already been established. Be- 
cause the fluxes only involve modes at comparable scales, 
they cannot transfer energy from very distant scales into 
scale £ within an inertial-inductive range. This is not 
true in the Batchelor range since the velocity field is very 
smooth there (<r u = 1), violating the condition ^ for 
scale-locality of energy flux in the infrared. 

The studies [H, 0, Q considered more traditional spec- 
tral transfers such as T ub (K,P) = (t>[ p J (b- V)u^) and 
T bu {K,P) = ( U [ p ](b-V)b[ K l), where and 
are spectrally band-passed fields for some interval of 
wavenumbers around K. Since T ub (K,P) = — T(, U (P, K), 
these can be interpreted (with some caution) as energy 
transfer from the velocity field in band [K] to the mag- 
netic field in band [P]. Is it possible for the dominant 
transfers to be between distant bands in an inertial- 
inductive range? The answer is no, if [K] is a dyadic 
(octave) wavenumber band [K/2,K]. It is necessary to 
use such bands, of equal width on a logarithmic scale, 
in order to permit simultaneous localization of modes in 
Fourier and physical space (within the limits of the un- 
certainty principle) . We note that this is crucial for phe- 
nomenological arguments based upon Alfvenic wavepack- 
ets with both size and wavenumber specified. The condi- 
tions which replace ^ are, for a 



u, b with < o-p < 1 



(|a 



[K]w\l/p 



K' 



Va 



[K]\P\l/p 



(7) 



bee If P < K/2, then wavenumber conservation 

implies that T ub (K,P) = -(u^(b^ K / 2 ~ p ' K + p ^V)h^}. 
Using the Holder inequality, this expression is bounded 

by (|Vb[ P ]| 3 ) 1 /3(| u [^|3)l/3 ( | b [K/2- W P]|3 ) l/3_ By 



(Tn (Tn 



|T„ 6 (^,-P)| < (constOP 1 -^ 

Since o\ < 1, such transfers for P <C K are negligible. 
For P > 2K, T ub (K,P) = (b^(bW 2 - K ' P + K ^-V)u^), 
so Holder inequality and ([7]) imply 



\T ub {K,P)\ < (const.) K 1 -^ P 



-2<7% 



Since 03 > 0, transfers for P ^> K are also negligible. 

To test these conclusions numerically we analyze a time 
snapshot of our 1024 3 MHD simulation in the statistical 
steady state. The kinetic and magnetic energy spectra 
of the flow, have a reasonable power-law scaling until 
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FIG. 



The transfers | (djuf'Bf 1 Bj) | (o) and 
(x). Straight lines have ±2/3-slopes and 



extend over the fitting range, which yields a decay rate of 



for (o) and 



for (x). Inset shows velocity 



(solid line) and magnetic (dashed-dotted line) energy spectra, 



which scale close to E n 



E b ~ k~ 



over k £ [5, 




[P/2.P] R [100,200] 



Bj) (solid line, 
Bj) 



[P-l,P]^j[199,200] 



FIG. 2: The transfers (djuf / *'*' s B i 
same as (o) plot in Fig.[TJ and 10 3 x (dj 
(dashed line). The latter is multiplied by 1000 for comparison. 

around k = 80 (inset to Fig.QJ. The transfers plotted 
in Fig. [I] exhibit off-diagonal (P ^ K) decay close to our 
rigorous upper bounds with exponent a\ = 1/3. However, 
the value of this exponent determined from our numerical 
data (not shown) is closer to erg = 1/4, consistent with 
the predictions of P, 0, 0]. For this value we obtain 
rigorous upper bounds O(P - 75 ) for P < K and 0(P~ 5 ) 
for P 3> K, which are also close to the observed scaling. 

How are our exact results to be reconciled with the 
recent numerical studies that reach the opposite conclu- 
sion? A full discussion is given in our longer work [lj|, 
but we make a few remarks here. 0, Q discussed sim- 
ulations at lower resolution than ours without carrying 
out a systematic scaling analysis. As for 0], they had an 
anomalously strong strain at the forcing scale L, which 
can dominate over the local strain at scales i < L in an 
inertial- inductive range of limited extent. We also ob- 
serve this effect over a finite range if we permit such an 
"energy spike" at the forcing scale, but it becomes weaker 
as the amplitude of the spike decreases or as the length 
of the power-law scaling range increases. Finally, [|| ap- 



pealed to spectral transfers to justify their claim that the 
magnetic field at scales I in the inertial-inductive range 
receives energy from straining motions at all larger scales 
> £, especially from scale L ^> I. However, their DNS 
study used Fourier bands of linear size [K — 1,K], which 
correspond to plane-wave modes which are nonlocalized 
in space, unlike the Alfven wavepackets employed in phc- 
nomenological arguments. Such bands do not properly 
account for the exponentially growing number of local 
triads at higher wavenumbers, whose aggregate contribu- 
tion dominates transfers defined with logarithmic bands 
[l3j |. Fig.[2]reproduces the numerical result of Fig. 8 of @ 
(dashed-line), together with our own DNS results using 
log-bands. Clearly, the nonlocal effects observed by Q 
represent miniscule amounts of energy transfer compared 
with the net contribution of local triads and become even 
smaller as the scale range increases. In short, the numer- 
ical results in 0, H, 0, H| do not support any asymptotic 
non-locality of energy cascade in MHD turbulence. 
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